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A carbon nanotube is an ideal object for understanding 
the atomic scale aspects of interface interaction and friction. 
Using molecular statics and dynamics methods different types 
of motion of nanotubes on a graphite surface are investigated. 
We found that each nanotube has unique equilibrium orien- 
tations with sharp potential energy minima. This leads to 
atomic scale locking of the nanotube. The effective contact 
area and the total interaction energy scale with the square 
root of the radius. Sliding and rolling of nanotubes have dif- 
ferent characters. The potential energy barriers for sliding 
nanotubes are higher than that for perfect rolling. When 
the nanotube is pushed, we observe a combination of atomic 
scale spinning and sliding motion. The result is rolling with 
the friction force comparable to sliding. 



Although the fundamental aspects of friction have been 
studied for more then centuries, our knowledge about its 
microscopic aspects is very limited |l|. The invention 
of atomic force microscope (AFM) and its applica- 
tion in measurements of atomic scale friction |3|] (friction 
force microscope - FFM) have made a great impact on 
the studies of friction. A carbon nanotube is a stable 
nanoobject having cylindrical shape ||, thus ideal for 
understanding atomic scale friction. M. Falvo et. al. 
showed that it is possible to slide, rotate and roll car- 
bon nanotubes on a graphite surface M. They demon- 
strated that a nanotube has preferred orientations on the 
graphite surface and prefer rolling than sliding when it is 
in atomic scale registry with the surface. In this study, we 
carried out molecular statics, dynamics calculations and 
studies of stick-slip motion for a variety of nanotubes. We 
found that: (i) A nanotube has sharp potential energy 
minima leading to orientational locking. The locking an- 
gles are directly related to the chiral angle, (ii) Sliding 
and rolling of nanotubes have different characters. The 
energy barriers for sliding are higher then the barriers for 
perfect rolling, (iii) The effective contact area and total 
interaction energy scale with the square root of the ra- 
dius, (iv) A combination of sliding and spinning motion 
is observed when the tube is pushed. The net result is 
rolling with the friction force comparable to the corre- 
sponding force for sliding. 

The character of interaction between the moving object ( 
atom, molecule or any nanoparticle ) and the underlying 
surface defines the motion ||. The interaction energy 
may consist of short-range, attractive interaction energy 



due to chemical bonding; short-range repulsive energy 
and long-range, attractive van der Waals energy. The 
interaction between a carbon nanotube and a graphite 
surface is similar to that between two graphite planes 
which is weak and van der Waals in origin. To investi- 
gate the overall behavior of the motion of a carbon nan- 
otube on a graphite surface, we represent the interaction 
between the tube and the graphite surface atoms by an 
empirical potential of Lennard- Jones type |7]] which was 
used extensively to study solid Cqq and nanotube J3|. 
Recent theoretical calculations [[| showed that multiwall 
nanotubes on a graphite surface are not deformed signif- 
icantly. The atomic scale motion is determined mostly 
by the interaction of the outmost layer of the nanotube 
with the surface. In this work we studied rigid single wall 
nanotubes with different chiralities and radii. 
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FIG. 1. The interaction energy as a function of rotation 
angle between the nanotube axis and the graphite lattice for 
(10,10), (30,0) and (20,10) nanotubes. Each nanotube has 
unique minimum energy orientations repeating in every 60 
degrees. 

The energy barriers related to the motion of nanotubes 
can be conveniently analyzed by calculating the variation 
of the potential energy, Ep, and corresponding force dur- 
ing the motion. Four different types of motion are consid- 
ered: spinning, rotating, sliding and rolling. During each 
step of motion the height of the tube is optimized. We 
first spin and rotate the nanotubes in order to find the 
equilibrium positions. Fig. 1 shows the interaction en- 
ergy as a function of the rotation angle between the tube 



axis and the graphite lattice (All the data given in this 
study is for per A length of the nanotubes) . Each nan- 
otube has unique equilibrium orientations repeating in 
every 60°, reflecting the lattice symmetry of the graphite. 
The variation of energy near the minimum is very sharp 
which causes atomic scale locking of the nanotube. Lock- 
ing angles are different for different nanotube and they 
are the direct measure of the chiral angle. This provides 
a novel method for measuring the chiralities of carbon 
nanotubes. Other important point in Fig. 1 is that the 
energy variation between two consecutive energy minima 
is very small (except near the minima) . Thus the force 
needed to rotate a nanotube is very small when the tube 
is out-of-registry. These results are in good agreement 
with the recent experiments by M. Falvo et. al. 
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FIG. 2. (a) The variation of the interaction energy,_Bp of a 
(20,10) nanotube as a function of sliding and rolling distances, 
(b) The friction force in stick-slip motion of the (20,10) nan- 
otube when attached to an AFM tip with spring constant 
~ 8.0 x 10~ 3 N/m. 

Next, we studied sliding and rolling of carbon nan- 
otubes. In Fig. 2(a) the variation of the total interac- 
tion energy, Ep(s) of a (20,10) nanotube as a function 
of sliding distance, s is shown. The energy variation is 
very small for rolling since perfect atomic scale registry 
is always maintained in the contact region. On the other 
hand in sliding motion, all the atoms in the contact re- 
gion move simultaneously out-of-registry to higher energy 
positions and contribute to the energy barrier of motion. 
When the nanotube is attached to an AFM tip, stick- 
slip motion occurs. The tube first sticks and then slips 
suddenly (slides or rolls or both) when the force exerted 
by the tip is sufficiently large [Q. The friction force in 



stick-slip motion of the nanotube when it is attached to 
an AFM tip ( with string constant ~ 8.0 x 10~ 3 N/m ) is 
shown in Fig. 2(b). The area in the hysteresis curve gives 
us the amount of energy dissipated during the stick-slip 
motion. 

To investigate the tube dependence, we performed cal- 
culations with tubes in different chiralities or radii. Fig. 
3(a) shows the interaction energy as a function of the 
nanotube radius, R. We found that the effective contact 
area and the interaction energy scale with the square 
root of the radius of the nanotube. The interaction en- 
ergy is independent of chirality. In-registry sliding force 
(when the tube is in a minimum energy orientation) is 
also scale with the square root of the radius. However, 
the sliding force is different for nanotubes with different 
chiralities and the same radii. Fig. 3(b) shows the force 
for in-registry sliding and rolling. For a typical nanotube 
( radius ~ 13 nm, length ~ 600 nm ), the sliding force 
value is estimated as ~ 87 nN for an armchair tube and 
~ 43 nN for a zigzag tube in good agreement with the 
friction force values measured in the experiments B ( 
- 50 nN ). 
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FIG. 3. (a) The interaction energy as a function of square 
root of the nanotube radius, \> r R. The filled circles correspond 
to armchair tubes and the hollow circles correspond to zigzag 
tubes with different radii, (b) The corresponding force for 
in-registry sliding and rolling of the nanotubes as a function 
ofVJ?. 

The static calculations we discussed gives insight of ideal 
sliding and rolling. However, the dynamical behavior is 
crucial for the competition between sliding and rolling in 
the course of motion. In our model, the graphite sub- 
strate and the nanotube are considered to be rigid but 



the nanotube as a whole is able to move having transla- 
tional and spinning degrees of freedom. Constant lateral 
force was applied to the nanotube for a short period of 
time ( 50 ps ) and then the motion of the nanotube was 
analyzed. In the same way, we applied constant torque 
or combinations of torque and lateral force to the nan- 
otube. Total energy of the system was kept constant. 




FIG. 4. (a) After initial push, sliding and spinning compo- 
nents of the kinetic energy of the (10,10) nanotube as func- 
tions of time are represented by the solid and dashed lines 
respectively. Notice the switching between spinning and slid- 
ing motions in the atomic scale, (b) The total force acting on 
the tube in the direction of motion (c) The sliding and spin- 
ning distances (angle multiplied by R) as functions of time. 

After constant lateral force is applied on the nanotube, 
slide-spin motions in the atomic scale are observed. 
When the atoms in the contact region are in atomic scale 
registry it is easier for the nanotube to slide. Then the 
nanotube atoms move from in-registry to out-of-registry 
positions and it is easier for the nanotube to spin. By 
spinning the nanotube decreases its potential energy and 
the atoms recover the in-registry positions. The switch- 
ing of the tube motion between spinning and sliding can 
be clearly seen in the sliding and the spinning component 
of the kinetic energy shown in Fig. 4(a). The switching 
between spinning and sliding is in the atomic scale and 
directly related to the corrugation of interaction energy. 
The total force acting on the tube in the direction of 
motion is shown in Fig 4(b). The maximum value of the 
force is comparable to the force required to slide the nan- 
otube (see Fig. 3(b)). Sliding and spinning distances ( 



angle multiplied by R) as functions of time are shown 
in Fig. 4(c). Ideal rolling would be a perfect overlap of 
sliding and spinning distances meanwhile slide-spin mo- 
tion gives oscillations and the net result is equivalent to 
rolling. 

For a better understanding of the slide-spin motion, we 
plot the interaction energy as a function of sliding dis- 
tance and spinning angle in Fig. 5(a). The trajectory 
for ideal rolling is a line at the bottom of the valley like 
regions. However, a nanotube performing slide-spin mo- 
tion follows an oscillating path in these valleys (see Fig. 
5(b)) with an amplitude of oscillation depending on the 
initial kinetic enerev. 
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FIG. 5. (a) The variation of the interaction energy as 
a function of sliding distance and spinning angle, (b) The 
trajectories correspond to the slide-spin motion of a nanotube 
having lower and higher total kinetic energies on the plane 
defined in (a) 

When the system is coupled to a heat bath or energy dis- 
sipation due to friction is considered, there are changes in 
the slide-spin motion. We modeled the dissipation by ad- 
ditional velocity dependent forces O] on the atoms close 
to the contact region. The results are presented in Fig. 6. 
Without energy dissipation the tube oscillates in the val- 
ley like regions of potential energy surface (see Fig. 5(a)) 
in the slide-spin motion. When the energy dissipation is 
considered the total kinetic energy decreased and there 
is more mixing between sliding and spinning. Eventu- 
ally, the nanotube performs ideal rolling. If the nanotube 
has very high kinetic energy, it slides over many surface 
unit cells. But due to energy dissipation the tube's to- 



tal kinetic energy decreases and slide-spin motion starts. 
Afterwards the motion is close to ideal rolling (As seen 
in Fig. 6(b)). This atomic scale picture of rolling is 
very similar to the rolling of macroscopic objects. Recent 
molecular dynamics simulations JL2| | with a full relaxed 
system by J. D. Schall and D. W. Brenner find similar 
conclusions. 
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FIG. 6. (a)Sliding and spinning component of the kinetic 
energy of the (10,10) nanotube as functions of time are rep- 
resented by the solid and dashed lines respectively (b) The 
nanotube's trajectories on the surface defined in fig. 5(a) 
for different initial kinetic energies. The trajectories from 
higher to lower initial kinetic energies are represented by dot- 
ted-dashed, solid and dashed lines, respectively. The kinetic 
energy values in (a) corresponds to the trajectory plotted by 
solid line. 

To conclude, we investigated different types of motion of 
carbon nanotubes on a graphite surface. Each nanotube 
has unique minimum energy orientations with respect to 
the surface structure. The variation of interaction en- 
ergy is very sharp leading to orientational locking of the 
nanotubes. The locking angles are direct measure of the 
chiral angles and this provides a novel method for mea- 
suring the nanotube chirality. We found that the effective 
contact area and the total interaction potential energy 
scale with square root of the radius of the nanotube. A 
combination of atomic scale spinning and sliding motion 
is observed when the nanotube is pushed. 
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